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The  indirect  band  gap  of  icosahedral  B12As2  (IBA)  has  been  determined  by  variable  temperature 
photoluminescence  measurements  (8K-294K)  on  solution-grown  bulk  samples.  In  addition, 
evidence  of  three  shallow  acceptor  levels  and  one  shallow  donor  level  is  reported.  The  low- 
temperature  spectra  were  characterized  by  broad  and  intense  deep  defect  emission,  donor-acceptor 
pair  (DAP)  bands,  and  exciton  recombination.  The  appearance  of  DAP  emission  verifies  the 
incorporation  of  a  donor  in  IBA,  which  has  not  been  reported  previously.  The  temperature 
dependence  of  the  free  exciton  (FE)  intensity  reflected  a  FE  binding  energy  of  45  meV.  The 
variation  of  the  FE  peak  position  with  temperature  was  fitted  with  both  Varshni  and  Passler  models 
to  determine  an  expression  for  the  temperature  dependence  of  the  indirect  band  gap.  The  resulting 
low  and  room  temperature  band  gaps  are  Eg(0)  =  3.470 eV  and  Eg(294  K)  =  3.373  eV,  respectively. 

The  latter  is  not  consistent  with  previous  reports  of  the  room  temperature  band  gap,  3.20  eV  and 
3.47  eV,  derived  from  band  structure  calculations  and  optical  absorption,  respectively.  The  origin 
of  these  discrepancies  is  discussed.  The  DAP  spectra  reveal  three  relatively  shallow  acceptors  with 
binding  energies  of  « 175,  255,  and  291  meV,  and  a  shallow  donor  with  binding  energy  «25  meV. 
Although  the  identity  of  the  individual  acceptors  is  not  known,  they  appear  to  be  associated  with 
the  light-hole  band.  The  small  donor  binding  energy  is  suggestive  of  an  interstitial  donor  impurity, 
which  is  suspected  to  be  Ni.  ©2012  American  Institute  of  Physics. 

[http://dx.doi.org/10.1063/L4729920] 


I.  INTRODUCTION 

Icosahedral  boron  arsenide,  B12As2  (IBA),  is  a  wide 
bandgap  semiconductor  that  exhibits  unusual  properties  as 
the  result  of  the  bonding  within  the  12-atom  boron  icosahe¬ 
dron  that  is  an  integral  part  of  its  crystal  structure.1"  One  of 
its  most  interesting  properties  is  the  ability  to  “self-heal” 
under  intense  ionizing  radiation  that  would  render  other  mate¬ 
rials  amorphous.2'1  Emin  has  argued"  that  this  radiation  hard¬ 
ness  results  from  the  unique  three-center  bonds  between 
icosahedral  boron  atoms.  The  bond  charge  is  retained  on  the 
icosahedron  when  a  boron  atom  is  displaced  due  to  the  radia¬ 
tion;  thus  resulting  in  a  Coulomb  attraction  between  the  va¬ 
cancy  and  interstitial  that  enhances  defect  recombination. 
IBA  is  also  a  hard  material  with  a  high  melting  temperature 
(>2000  °C).  This  unique  combination  of  physical  characteris¬ 
tics  makes  IBA  a  natural  choice  for  applications  that  must 
withstand  high  temperatures  or  intense  radiation,  such  as  ther¬ 
moelectrics,  neutron  detectors,  and  alpha-  and  beta-voltaics.2 

While  there  has  been  sparse  but  continuing  interest  in 
the  properties  of  IBA  over  the  years,  the  material  remains 
poorly  characterized  and  will  require  significant  progress  in 
materials  growth  to  realize  its  potential.  Small  bulk  crystals 
prepared  from  a  liquid  Pd-B  melt  were  grown  in  early  meas¬ 
urements  of  the  thermal  conductivity  of  IBA.4  Within  the 
last  several  years,  there  has  been  renewed  activity,  mostly 
focused  on  growing  epitaxial  layers  of  IBA  on  Si  (Ref.  5) 
and  on  4H-,  6H-  and  15R-SiC.6  8  Whiteley  et  al.9  have 
recently  investigated  the  growth  of  bulk  crystals  from  a  Ni-B 
solution.  In  all  of  these  studies,  the  focus  of  the  material 


characterization  has  been  on  structural  properties,  with  x-ray 
diffraction,7  transmission  electron  microscopy,7'810  and 
Raman  scattering1  1,2  as  the  primary  experimental  probes. 

Recently  Gong  et  al. 1 1  reported  on  the  fabrication  and  per¬ 
formance  of  an  IB  A/n-type  4H-SiC  p-n  junction  device,  where 
the  IBA  was  p-type  with  a  carrier  concentration  of  roughly 
1017cm-3.  In  fact,  as-grown  boron-rich  semiconductors  are 
generally  p-type.  This  occurs  because  inter-icosahedral  bond¬ 
ing  is  maintained  when  an  As  atom  is  replaced  by  a  vacancy  or 
either  a  B,  group  IV,  or  group  VI  atom,  and  the  resulting  hole 
on  the  As  site  leads  to  p-type  conductivity.2  For  example,  Xu 
et  a/.14  demonstrated  that  the  conductivity  of  high-resistivity 
p-type  IBA  epilayers  grown  on  semi-insulating  6H-SiC  could 
be  effectively  modulated  by  Si  doping,  resulting  in  p-type  con¬ 
duction  with  a  x  105  reduction  in  the  resistivity.  Hall  effect 
measurements  on  similar  as-grown  epilayers  («3  /(in  thick) 
also  exhibited  p-type  conductivity  with  hole  concentrations15 
in  the  range  («1013-1016)  cm-3.  Thus  far,  n-type  conductivity 
has  not  been  reported  in  IBA. 

Other  studies  of  IBA  have  focused  on  thermal  properties 
(Seebeck  coefficient,  thermal  conductivity16),  calculations 
of  the  vibrational  properties,17'18  and  band  structure 
calculations, 19-22  which  indicate  that  the  band  gap  of  IBA  is 
indirect,  with  a  conduction  band  minimum  at  the  A-point  of 
the  Brillouin  zone.  Optical  studies  have  been  limited  to 
Raman  scattering;23  early  work  by  Slack  et  al ,24  on  refrac¬ 
tive  index  and  optical  absorption,  where  the  indirect  band 
gap  was  reported  to  be  3.47  eV  at  room  temperature;  infrared 
(lattice)  absorption  appearing  in  Ref.  17;  and  recent  room 
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temperature  spectroscopic  ellipsometry  measurements  by 
Bakalova  et  a/.,21  which  provided  information  on  the  dielec¬ 
tric  function  and  interband  critical  points.  The  combination 
of  this  ellipsometry  data  with  first  principles  band  structure 
calculations21  enabled  critical  point  assignments  and  the 
determination  of  the  required  self-energy  correction 
(«0.6eV)  for  the  calculated  direct  band  gaps.  No  experimen¬ 
tal  data  were  available  near  the  indirect  gap  due  to  compet¬ 
ing  impurity  absorption.  Consequently,  the  fundamental 
band  edge  was  determined  by  adding  the  self-energy  correc¬ 
tion  obtained  for  direct  transitions  to  the  calculated  indirect 
gap  (2.56  eV),  leading  to  a  band  gap  of  «3.2eV.  This  band 
structure  calculation  also  provided  the  first  determination  of 
the  electron  and  hole  effective  masses  for  1BA.21 

Thus,  there  are  two  very  different  values  in  the  literature 
for  the  indirect  band  gap  of  IBA:  3.20  eV  (Ref.  21)  and 
3.47  eV.24  In  this  work,  we  address  this  issue  by  carrying  out 
photoluminescence  (PL)  measurements  over  a  broad  temper¬ 
ature  range.  It  is  somewhat  surprising  that  PL  spectroscopy 
has  not  been  reported  previously  for  IBA,  since  the  material 
has  been  of  interest  for  at  least  40  years.  Optical  characteri¬ 
zation  such  as  PL  can  provide  information  on  the  semicon¬ 
ductor  band  gap  and  on  the  defects  and  impurities  that  define 
or  limit  its  transport  properties.  Since  the  existing  values  for 
the  indirect  band  gap  were  measured  at  room  temperature, 
while  PL  spectral  features  are  more  easily  identified  at  low 
temperatures,  the  present  measurements  are  carried  out  over 
the  full  temperature  range  between  8  K  and  294  K. 

II.  EXPERIMENTAL 

The  IBA  crystals  were  grown  by  the  procedure  described 
in  detail  by  Whiteley  et  al.9  In  a  sealed  quartz  tube,  boron 
and  nickel  were  placed  in  a  boron  nitride  boat  at  one  end,  and 
arsenic  was  placed  at  the  other  end.  The  base  pressure  in  the 
quartz  tube  before  sealing  is  typically  5  x  10-6  torr.  To  begin 
the  crystal  growth,  the  zone  containing  boron  and  nickel  alloy 
was  heated  to  1 150  °C  within  1  h  and  held  at  this  temperature 
for  120  hours.  The  opposite  end  of  the  quartz  tube  that  origi¬ 
nally  contained  the  arsenic  was  slowly  heated  up  to  600  °C, 
where  it  was  held  for  120  hours.  Crystal  growth  was  initiated 
by  slowly  cooling  to  1000°C  at  a  rate  of  2°C/h.  The  furnace 
was  then  turned  off,  and  the  quartz  tube  was  allowed  to  natu¬ 
rally  cool  to  room  temperature. 

Material  growth  was  carried  out  with  both  low  and  high 
purity  source  materials.  Crystals  from  the  low-purity  growth 
generally  contained  significant  concentrations  of  metallic 
and  nonmetallic  impurities.  To  grow  the  high  purity  IBA 
crystals,  care  was  taken  to  use  all  high  purity  source  materi¬ 
als:  The  boron  was  99.9999%  pure,  isotopically  enriched 
boron- 11,  the  arsenic  was  99.99999%  pure  (metals  basis), 
and  the  high  purity  nickel  used  had  C,  O,  and  Si  concentra¬ 
tions  of  less  than  50,  70,  and  1  ppm  respectively. 

After  the  growth  process,  the  crystallized  B12As2  was  then 
removed  from  the  matrix  by  etching  away  the  nickel  in  aqua 
regia  at  room  temperature.  The  low  purity  IBA  crystals  were 
generally  dark  gray  or  black,  and  opaque.  Many  of  the  high  pu¬ 
rity  IBA  crystals  were  also  dark  gray,  or  black,  and  opaque, 
but  some  of  the  thinnest  crystals  were  clear  and  transparent. 


TABLE  I.  Impurity  concentrations  measured  in  IBA  crystals  by  glow  dis¬ 
charge  mass  spectrometry. 


Impurity 

Growth  from  low  purity 
sources:  Sample  C 

Growth  from  high  purity  sources: 
Samples  A  and  B 

Concentration 

Concentration 

ppm 

cm  3 

ppm 

cm~3 

Si 

400 

3.0  x  1019 

11 

8.3  x  1017 

Ni 

260 

9.4  x  10IS 

2200 

8.0  x  1019 

S 

89 

5.9  x  101S 

<0.5 

<3.3  x  1016 

Er 

25 

3.2  x  1017 

<0.15 

<1.3  x  1015 

Se 

12 

3.2  x  1017 

<1 

<2.7  x  1016 

P 

6.5 

4.5  x  1017 

1.4 

9.6  x  1016 

Three  samples  were  employed  in  this  study:  Samples  A 
and  B  were  grown  using  the  high  purity  sources  while  sam¬ 
ple  C  was  grown  using  the  lower-purity  sources.  The  impu¬ 
rity  levels  for  both  low  purity  and  high  purity  IBA  crystals 
were  measured  by  glow  discharge  mass  spectrometry  and  are 
shown  in  Table  I.  While  the  concentrations  of  most  impur¬ 
ities  dropped  dramatically  when  the  high-purity  sources 
were  used,  this  was  accompanied  by  about  an  order  of  mag¬ 
nitude  increase  in  the  incorporation  of  Ni. 

Photoluminescence  measurements  were  carried  out  with 
the  samples  placed  inside  a  variable  temperature  optical  cryo¬ 
stat.  Excitation  was  provided  by  <20  mW  of  325  nm  light 
from  a  He-Cd  laser.  The  emitted  radiation  was  collected  and 
focused  into  a  0.85-m  double  spectrometer  employing  1800 
gpm  gratings  and  was  detected  by  a  cooled  GaAs  photomulti¬ 
plier.  The  spectral  bandpass  was  usually  set  at  «2.2meV. 

III.  RESULTS  AND  DISCUSSION 

A.  Overview — low  temperature  photoluminescence 
spectra 

Low  temperature  PL  spectra  taken  over  a  broad  energy 
range  are  shown  in  Fig.  1  for  the  two  high-purity  samples  A 
(dotted  line)  and  B  (dashed  line)  as  well  as  for  the  low-purity 
sample  C  (solid  line).  The  dominant  spectral  features  are 


FIG.  1 .  Low  temperature  PL  spectra  of  the  two  high  purity  samples  A  (dot¬ 
ted  line)  and  B  (dashed  line)  and  the  low  purity  sample  C  (solid  line). 
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remarkably  similar  to  those  observed  in  other  wide  bandgap 
semiconductors,  such  as  GaN,  which  is  usually  dominated 
by  a  broad  and  intense  deep  defect  band,  a  donor-acceptor 
pair  (DAP)  region,  and  a  near-gap  region  dominated  by  exci- 
tonic  recombination.  In  the  IB  A  spectra  of  Fig.  1,  a  broad, 
deep  defect-related  band  is  observed  centered  at  2.4  eV.  In 
the  following  sections,  we  will  discuss  in  detail  the  excitonic 
recombination  from  IBA,  represented  in  Fig.  1  by  the  two 
sharp  emissions  at  3.309  and  3.366  eV  and  DAP  recombina¬ 
tion  observed  in  the  range  3. 2-3.4  eV  in  the  figure.  Compar¬ 
ing  the  spectra  of  the  two  high-purity  samples,  sample  A 
exhibits  much  stronger  defect  and  DAP  emission  and  much 
weaker  exciton  recombination  than  sample  B,  indicating  that 
sample  B  is  of  higher  quality.  The  low-purity  sample  C 
exhibits  only  a  very  intense  defect  band  and  no  exciton- 
related  emission,  which  is  consistent  with  a  very  high  con¬ 
centration  of  defects  and/or  impurities. 

The  PL  spectra  from  all  three  samples  exhibited  signifi¬ 
cant  spatial  inhomogeneity,  as  the  spectra  varied  substantially 
when  the  exciting  laser  spot  was  moved  over  the  samples. 
This  inhomogeneity  was  also  confirmed  through  PL  imaging 
measurements  and  was  evident  from  the  overall  spatial  varia¬ 
tion  in  the  spectra.  The  data  presented  in  the  rest  of  this  work 
were  always  obtained  from  the  high-quality  region  of  each 
sample,  characterized  by  a  maximum  exciton  emission  inten¬ 
sity  and  a  minimum  defect  band  intensity. 

Broad  and  intense  defect/impurity  bands  emitting  well 
below  the  band  gap  appear  to  be  a  common  feature  in  the  PL 
spectra  of  wide  band  gap  semiconductors.  The  width  and 
general  lineshape  of  the  «2.4eV  emission  observed  from 
IBA  is  reminiscent  of  the  «2.25  eV  “yellow  band”  in  GaN, 
which  is  thought  to  be  a  carbon-related  defect,25’26  and  the 
broad  «2.4eV  emission  in  n-type  4H-SiC  related  to  a  deep 
boron-complex  center.-7  Determining  the  nature  of  the 
defect  responsible  for  these  emissions  requires  a  study  of  the 
effect  of  impurity  and/or  defect  concentrations  on  the  PL 
spectrum.  While  these  measurements  have  not  been  carried 
out  in  IBA,  from  Table  I,  it  is  evident  that  sample  C  contains 
an  order  of  magnitude  less  Ni  and  considerably  higher  con¬ 
centrations  of  all  other  impurities  as  compared  to  samples  A 
and  B.  While  this  suggests  that  the  high  concentration  of  Ni 
is  not  responsible  for  the  quenching  of  the  exciton  emission 
in  sample  C,  it  is  more  likely  that  the  strong  exciton  PL  in 
samples  A  and  B  arises  from  the  high-quality  region  of  each 
sample,  which  contains  a  much  lower  Ni  concentration. 
There  is  also  a  large  variation  in  the  defect-  and  exciton- 
related  PL  intensities  of  the  two  “high  purity”  samples  A  and 
B,  which  makes  it  difficult  to  make  more  meaningful  correla¬ 
tions  between  the  relative  PL  intensities  and  the  level  of 
impurities  in  these  samples. 

B.  Bound  exciton  spectra 

At  very  low  temperatures,  optical  excitation  produces 
photoexcited  electron-hole  pairs  that  rapidly  form  free  exci- 
tons  which  are  just  as  rapidly  captured  at  (usually  neutral) 
impurities,  forming  bound  excitons  (BEs).  Until  the  tempera¬ 
ture  is  raised,  the  BE  state  is  the  most  stable  configuration 
for  the  carriers.  The  bound  exciton  can  annihilate  by  recom¬ 


bining  radiatively,  as  BE  PL,  or  it  can  dissociate  thermally, 
leaving  a  neutral  impurity  and  either  a  free  exciton  or  a  free 
electron-hole  pair.  While  these  two  dissociation  mechanisms 
are  most  common,  others  involving  either  Auger  recombina¬ 
tion  and/or  a  change  in  the  charge  state  of  the  impurity  are 
also  possible.28-30  As  we  do  not  have  reliable  information  of 
the  carrier  type  or  compensation  level  of  the  samples  investi¬ 
gated  here,  the  bound  excitons  can  be  associated  with  neutral 
or  possibly  ionized  donors  or  acceptors. 

The  dominant  sharp  line  at  3.366  eV  in  Fig.  1  is  identi¬ 
fied  as  an  exciton  bound  to  a  shallow  impurity,  or  “shallow 
bound  exciton”  (SBE),  based  not  only  on  its  relatively  nar¬ 
row  line  width  and  its  position  in  the  spectrum  but  also  from 
its  temperature  dependence,  which  is  discussed  in  Sec.  Ill  C. 
The  SBE  line  is  shown  in  more  detail  as  the  solid  line  in  Fig. 
2(a).  Because  IBA  is  an  indirect  band  gap  semiconductor,  it 
was  not  initially  clear  whether  this  intense  emission  band 
was  a  zero-phonon  line  (ZPL)  or  resulted  from  a  phonon- 
assisted  transition.  However,  a  very  weak  emission  suspected 
of  being  the  ZPL  for  this  band  was  observed  at  3.419  eV, 
shown  magnified  by  500X  in  Fig.  2(a).  To  verify  its  origin, 
the  slow-varying  background  was  subtracted  from  the 


Energy  (eV) 


FIG.  2.  (a)  Low  temperature  PL  spectrum  of  the  3.366  eV  SBE  (solid  line) 
compared  to  the  scaled  and  shifted  weak  emission  at  3.419  eV  (dots),  (b) 
Log-intensity  plot  of  the  low  temperature  PL  in  the  exciton  region.  The 
DBE  exhibits  two  phonon  replicas  separated  by  Ep2  =  71.6meV,  while  a 
weak  SBE  ZPL  and  a  strong  1-P  replica  at  3.366  eV  are  observed. 


013508-4 


Klein  etal. 


J.  Appl.  Phys.  112,  013508  (2012) 


spectrum,  the  intensity  was  scaled  to  that  of  the  BE  line,  and 
the  band  was  shifted  down  in  energy  by  EPi  =  52.5  meV  to 
coincide  with  the  BE.  The  resulting  spectrum  (dots)  in  Fig. 
2(a)  is  almost  identical  to  that  of  the  BE  line  (solid  line),  ver¬ 
ifying  that  this  weak  emission  is  the  SBE  ZPL  and  that  the 
dominant  emission  is  a  one-phonon  (1-P)  replica  involving  a 

52.5  meV  («423cm-1)  momentum-conserving  (MC) 
phonon. 

Hence,  the  main  SBE  PL  transition  appears  to  be  phonon- 
assisted.  Dipole  selection  rules  require  that  the  MC  phonons 
be  of  odd  symmetry.  Fan18  has  recently  calculated  phonon  dis¬ 
persion  relations  for  IB  A.  Of  the  42  possible  modes,  10  infra¬ 
red  active  “u”  modes  were  identified:  4A2u  +  6EU.  One  of  the 
zone  center  Eu  modes,  calculated  at  486  cm- 1  and  observed  in 
infrared  absorption  at  483  cm”1,31  was  found  by  Fan  to  split 
into  two  modes  along  the  T-A  direction,  appearing  near 
506cm-1  and  406  cm-1  at  the  A-point.  Although  BE  recombi¬ 
nation  does  not  always  require  phonon  assistance,  in  this  case, 
it  appears  that  the  lower  of  these  two  A-point  phonons  may 
well  be  involved. 

In  Fig.  2(b),  the  spectrum  of  the  full  exciton  region  is 
shown  on  a  log-intensity  scale,  so  that  all  significant  features 
can  be  more  clearly  seen.  After  the  intense  SBE  line  at 
3.366  eV,  the  next  dominant  feature  is  located  at  3.309  eV, 
57  meV  lower  in  energy.  This  feature  can  be  seen  in  Fig. 
2(b)  to  have  two  further  phonon  replicas,  both  separated  by 
about  Ep2  =  71.6meV  («578cm-1).  If  the  3.309eV  band  is 
a  one-phonon  replica,  then  a  ZPL  would  be  expected  at 

71.6  meV  higher  energy,  at  3.381  eV.  This  occurs  near  the 
high  energy  edge  of  the  SBE  emission,  as  indicated  in  Fig. 
2(b).  After  careful  study  of  this  spectral  region,  no  PL  feature 
was  detected  at  this  energy,  suggesting  that  the  DBE  line  is  a 
ZPL.  As  we  will  show  in  Sec.  Ill  D,  this  corresponds  to  a 
DBE  about  116meV  below  the  exciton  gap.  If  the  transition 
were  phonon  assisted,  this  would  be  reduced  by  EP2  to 
«44meV.  These  are  large  BE  localization  energies.  Such 
deep  bound  excitons  (DBEs)  are  often  associated  with  transi¬ 
tion  metal  impurities.32-34  Thus,  the  high  Ni  content  in  these 
samples  suggests  the  possibility  that  the  DBE  emission  could 
result  from  excitons  bound  to  a  deep  Ni-related  center.  The 
DBE  line  also  exhibits  weaker  shoulders  at  lower  and  higher 
energy.  While  the  higher  energy  feature  is  close  to  one  pho¬ 
non  energy  below  the  SBE  line  and  could  be  a  second  pho¬ 
non  replica,  the  lower  feature  remains  unidentified. 

The  SBE  line  at  3.366eV  participates  with  a  423cm-1 
phonon  in  the  indirect  transition,  while  the  DBE  emission 
involves  a  578  cm- 1  phonon.  Because  the  two  BEs  experience 
very  different  carrier  localizations,  it  is  not  surprising  that 
they  couple  effectively  to  different  phonons.  While  an  Eu  lat¬ 
tice  mode,  calculated  by  Fan18  at  598  cm-1  and  observed  in 
absorption31  at  580  cm-1  could  be  associated  with  the  DBE 
indirect  transition,  it  is  also  possible  that  the  phonon  coupling 
to  the  DBE  is  a  local  vibrational  mode. 

The  optical  transitions  described  in  Fig.  2(b)  are  sum¬ 
marized  in  the  level  diagram  of  Fig.  3  as  the  wavy  lines,  with 
MC  phonons  indicated  by  the  saw-toothed  lines.  The  figure 
elucidates  the  relationship  between  the  observed  PL  transi¬ 
tion  energies  and  the  significant  energy  levels  of  the  system. 
The  DBE  emission  is  represented  as  a  ZPL. 


FIG.  3.  Schematic  level  diagram  indicating  optical  transitions  as  wavy  lines 
and  MC  phonons  as  saw-toothed  lines.  The  most  prominent  PL  features  in 
the  exciton  region  are  indicated,  as  well  as  low  temperature  PL  peak  posi¬ 
tions  and  the  characteristic  exciton  energy  separations  Afe,  Abe*  and  ADbe- 
The  DBE  line  at  3.309  eV  is  assumed  to  be  a  ZPL. 

C.  Temperature  dependence  of  the  exciton  spectra 

The  temperature  dependence  of  the  PL  spectra  in  the 
exciton  region  is  shown  in  Fig.  4  for  sample  B  over  the  tem¬ 
perature  range  8K-294K.  The  maximum  intensity  of  each 
spectrum  is  normalized  to  the  same  height  and  shifted  verti¬ 
cally,  so  that  relative  intensity  variations  can  be  compared 
directly.  At  the  lowest  temperatures,  the  spectra  exhibit  only 
the  DBE  and  SBE  emissions.  As  temperature  is  raised,  the 


FIG.  4.  Temperature  dependence  of  the  exciton  PL  spectra  of  sample  B.  All 
spectra  are  normalized  to  the  same  height  and  are  shifted  vertically. 
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FIG.  5.  Temperature  dependence  of  the  integrated  PL  intensities  of  the  (a) 
FE,  (b)  SBE,  and  (c)  DBE.  The  fits  of  the  data  to  Eqs.  (3),  (4),  and  (5), 
respectively,  appear  as  the  solid  lines.  The  inset  in  (b)  extends  to  a  broader 
temperature  range. 

absolute  intensity  of  the  SBE  decreases  slowly  while  that  of 
the  DBE  begins  to  increase  (see  e.g.  Fig.  5),  so  that  their  rela¬ 
tive  intensities  begin  to  change.  As  the  temperature  increases 
beyond  25  K,  a  new  band  emerges  from  the  high  energy  tail 
of  the  SBE  line  and  increases  in  intensity  relative  to  the  SBE 
as  temperature  is  further  increased,  becoming  totally  domi¬ 
nant  above  about  120K.  This  is  typical  behavior  for  FE  emis¬ 
sion:  At  very  low  temperatures,  the  BE  state  is  more  stable. 
Therefore,  unless  the  concentration  of  impurities  that  bind 
excitons  is  very  low,  FE  emission  at  very  low  temperature  is 
weak.  As  temperature  is  raised,  FEs  are  thermally  dissociated 


from  the  BEs,  thus  increasing  the  FE  concentration  and  the 
associated  FE  PL  intensity  and  decreasing  the  corresponding 
SBE  intensity.  At  higher  temperatures,  the  dissociation  of  the 
FE  into  free  carriers  dominates  and  the  FE  intensity  rapidly 
decreases  with  temperature.  This  dependence  of  the  FE  inten¬ 
sity  on  temperature  consequently  provides  a  measure  of  the 
FE  binding  energy,  AFE.  Once  AFE  is  determined,  the  meas¬ 
ured  temperature  dependence  of  the  energy  position  of  the  FE 
PL  peak  (1-P  replica),  lmFE_1P,  provides  the  information  nec¬ 
essary  to  determine  the  temperature  dependence  of  the  IBA 
bandgap  (see,  e.g.,  Fig.  3): 

Eg(T)  =  [1iz/fe_ip(T)  —  '/zkT]  +  AFE  +  EPi  =  Egx(T)  +  AFE, 

(1) 

where  EPJ  is  the  energy  of  the  MC  phonon  («52.5  meV/ 
423  cm-1).  Egx(T)  =  ImFE_1P(T)  —  '/zkT  +  EP1  is  the  exciton 
gap  energy  (i.e.,  the  position  where  the  ZPL  of  the  FE  emis¬ 
sion  would  appear).  The  factor  '/zkT  comes  from  the  FE  PL 
lineshape,  IFE(fm)  a  (hu  —  Egx)1/,2exp[— (hz/  —  Egx)/kT],  which 
peaks  '/ikT  higher  in  energy  than  the  actual  exciton  gap. 

The  SBE  is  localized  relative  to  the  FE  by  the  SBE 
localization  energy  ABE.  Since  the  peak  of  the  FE  PL  band  is 
at  '/ikT  higher  energy  than  the  exciton  gap,  the  separation  of 
the  FE  and  BE  PL  peaks  is  given  by 

[hj/FE_ip(T)  —  hz/BE-ip(T)]  =  Abe  +  V^kT.  (2) 

Consequently,  we  can  obtain  a  measure  of  ABE  from  the  sep¬ 
aration  of  these  two  peaks.  In  Sec.  Ill  E  we  find  that 
ABE  =  6.3  ±  1  meV,  independent  of  temperature. 

To  extract  the  temperature  dependence  of  the  DBE, 
SBE,  and  FE  integrated  intensities  and  peak  positions  from 
the  overlapping  spectral  features  in  Fig.  4,  the  data  were  fit¬ 
ted  using  commercial  peak  fitting  software. 

D.  Temperature-dependent  intensities 

The  temperature  dependent  free  and  bound  exciton  inten¬ 
sities  are  determined  by  the  steady-state  balance  of  photoex- 
cited  carrier  generation,  radiative  recombination,  and  the 
capture  and  thermal  liberation  of  excitons  from  their  bound 
states.  Ping35  has  worked  out  a  detailed  model  for  the  steady- 
state  case,  providing  temperature-dependent  expressions  for 
the  FE  and  SBE  populations  and  PL  intensities.  These  are  of 
the  form 

BE  (  )  ^  j  _|_  AbeAT  _|_  gg— (Abe+AfeIAT  ’ 

1  +  Ce-ABE//kT 

^  1  -)-  Ae~AfiE/kT  Be~(ABE+AFE)/kT  ^  ^ 

The  two  ionization  energies  appearing  in  the  exponents  in  Eqs. 
(3)  and  (4),  correspond  to  the  thermalization  of  the  BE  from 
the  impurity,  producing  either  a  FE  (ABE  term)  or  a  free 
electron-hole  pair  (ABE  +  AFE  term).  The  A’s,  B’s,  and  C’s 
involve  coefficients  for  carrier  emission  and  radiative  recom¬ 
bination  from  the  various  levels  of  the  system.  The  coupling  of 
all  three  levels  (free  carriers,  FEs  and  BEs)  through  emission 
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and  capture  is  responsible  for  the  appearance  of  both  exciton 
binding  energies,  ABE  and  AFE,  in  the  expression  for  the  inten¬ 
sity  of  each  exciton  emission.  The  measured  temperature  de¬ 
pendence  of  the  spectrally  integrated  PL  intensities  for  the  FE 
and  SBE  are  shown  as  the  filled  circles  in  Figs.  5(a)  and  5(b), 
respectively.  At  low  temperatures,  the  FE  intensity  initially 
increases  with  temperature  as  FEs  are  thermalized  from  BE’s, 
while  the  BE  intensity  correspondingly  decreases  monotoni- 
cally.  At  the  highest  temperatures,  both  excitons  are  thermal¬ 
ized  directly  into  free  electron-hole  pairs.  Both  sets  of  data 
(FEs  and  BEs)  can  be  fitted  (solid  lines)  with  Eqs.  (3)  and  (4) 
using  the  same  two  exciton  binding  energies:  ABE  =  6.3 
±  1  meV  and  AFE  =  45  ±  5  meV.  The  BE  localization  energy 
Abe  was  not  a  fitting  parameter  but  was  determined  directly 
from  the  separation  of  the  FE  and  BE  emission  peaks,  as  noted 
in  reference  to  Eq.  (2)  and  described  explicitly  in  Sec.  Ill  E. 
The  fitting  parameters  (beyond  an  overall  scaling  factor)  for 
the  SBE  were  A  =  4.5,  B  =  2600  and  for  the  FE  were 
A  =  0.005,  B  =  1400,  and  C  =  3000.  The  FE  binding  energy 
of  45  me V  is  somewhat  lower  than  the  66  meV  reported  by 
Bakalova  et  al ,,36  determined  from  the  effective  masses  eval¬ 
uated-1  from  first  principles  band  structure  calculations.  This 
is  discussed  in  detail  in  Sec.  Ill  F. 

The  temperature  dependence  of  the  DBE  intensity  is 
shown  in  Fig.  5(c).  The  PL  intensity  increases  with  increasing 
temperature  until  about  60  K  and  then  decreases,  qualitatively 
following  the  FE  population,  as  in  Fig.  5(a).  This  is  suggestive 
of  a  DBE  that  forms  by  trapping  FEs  that  have  been  thermal¬ 
ized  from  a  SBE  site.  The  DBE  can  also  sequentially  trap  free 
carriers.  We  have  extended  Ping’s  approach  slightly  to  add  the 
DBE  state  to  the  steady-state  model.  Experimentally,  this 
emission  peak  was  shifted  below  the  SBE  peak  by  ADBE 
=  h//SBF  —  hf/|)|iF  =  57  meV.  The  steady-state  DBE  popula¬ 
tion  is  determined  by  radiative  emission  and  by  transitions 
between  the  DBE  state  and  both  free  excitons  and  free  carriers. 
Incorporating  this  into  the  model  with  the  DBE  line  assumed 
to  be  a  ZPL,  the  temperature  dependence  of  the  DBE  popula¬ 
tion  (and/or  intensity)  can  be  approximated,  assuming  a  slow- 
varying  free  carrier  concentration,  in  the  form: 

Idbe  (T) 

1  +  Anpg(T) 

X  [  _|_  gg~(^BE+ADBE  +  Epi)/kT  _|_  Ce  fAiji  - A[jiji;  -  Hpi-  A|.E;)/k'l' 

(5) 

Here,  nFE(T)  represents  the  FE  population  (e.g.,  Eq.  (3)), 
while  the  two  transition  energies  appearing  in  the  exponents 
reflect  the  thermalization  of  the  DBE  state,  leading  (see,  e.g., 
Fig.  3)  to  the  liberation  of  either  FEs  (Abe  + Adbe  +  Ep1)  or 
free  electron-hole  pairs  (ABE  + Adbe  +  Ep1  +  AFE).  The 
expression  (Abe+ Adbe  +  Ep1)  corresponds  to  the  DBE 
localization  energy  Egx  —  ImBE,  «  116meV.  The  experimen¬ 
tal  data  (filled  circles)  are  shown  in  Fig.  5(c)  fitted  to  Eq.  (5) 
(solid  line)  using  the  same  values  of  ABE  =  6.3meV  and 
Afe  =  45  meV  as  determined  above  and  with  ADBE  =  57  meV 
and  EP1  =  52.5  meV  determined  directly  from  the  separation 
of  spectral  peak  positions.  The  fitting  parameters  were 
A=  1.54  x  10“3,  B  =  20,  and  C  =  3000. 


If  the  alternative  view  is  taken  that  the  3.309  eV  DBE 
PL  is  a  1-P  replica,  with  ZPL  at  3.381  eV  (not  observed), 
then  the  DBE  localization  energy  would  be  reduced  by 
EP2  =  71.6  me V  to  «44  meV — which  is  still  a  large  BE  bind¬ 
ing  energy.  As  the  temperature  dependence  in  Fig.  5(c)  is 
dominated  by  nFE(T),  a  similar  fit  of  the  data  may  be 
obtained  with  this  reduced  DBE  localization.  Consequently, 
the  lack  of  an  observed  ZPL  results  in  some  ambiguity 
regarding  the  DBE  localization  energy,  and  both  values  are 
entered  in  the  summary  of  IBA  material  parameters  appear¬ 
ing  in  Table  II. 

With  the  FE  binding  energy  of  45  ±  5  meV  measured 
from  the  temperature  dependent  FE  PL  intensity,  the  temper¬ 
ature  dependence  of  the  indirect  band  gap  can  now  be 
considered. 

E.  Temperature  dependence  of  the  energy  gap 

The  temperature  dependence  of  the  energy  positions  of 
the  FE,  SBE,  and  DBE  determined  from  the  peak-fitting 
analysis  of  the  temperature  dependent  spectra  is  shown  plot¬ 
ted  in  Fig.  6(a)  as  the  filled  triangles,  squares,  and  circles, 
respectively.  It  is  apparent  that  the  peak  position  of  the  exci¬ 
ton,  unlike  that  of  the  SBE  and  FE,  is  unaffected  by  the  tem¬ 
perature  dependence  of  the  band  gap,  up  to  temperatures 
greater  than  100  K.  This  is  behavior  consistent  with  a  exciton 
that  is  strongly  coupled  to  a  deep  binding  center:  Patrick  and 
Choyke7  investigated  DBEs  due  to  Ti  impurities  in  several 
SiC  polytypes  and  interpreted  the  independence  of  the  DBE 
recombination  energy  from  the  position  of  the  band  edge  as 
characteristic  of  excitons  that  involve  the  localized 
d-electrons  of  the  transition  metal.  Similarly,  Kimoto  et  al.34 
reported  Ti-related  DBE  emission  in  Ti:6H-SiC  that  exhib¬ 
ited  no  variation  in  spectral  position  over  a  wide  temperature 
range  ( 14-77  K).  Similar  behavior  would  be  expected  from 
excitons  bound  to  other  transition  metals.  Thus,  the  tempera¬ 
ture  independence  of  the  DBE  energy  suggests  the  possibil¬ 
ity  of  a  Ni-related  DBE. 

At  the  bottom  of  the  figure  the  measured  energy  differ¬ 
ences  between  the  FE  and  SBE  peak  energies  are  plotted  as 


TABLE  II.  Experimental  IBA  material  parameters. 


Parameter 

Description 

Value 

Eg(0) 

Low  temperature  gap 

3.470  ±  0.005  eV 

Eg(294  K) 

Room  temperature  gap 

3.373  ±  0.01  eV 

Eg(T) 

Temp.  dep. — Varshni 

3.470  ±  0.005  - 
4  x  10-4T2/(640  — T)  eV 

Temp.  dep. — Passler 

3.470  ±  0.005  - 

0.90{  [1  +  (T/450)]27  -  1  )(1/277)  eV 

Abe 

BE  localization  energy 

6.3  ±  1  meV 

Afe 

FE  localization  energy 

45  ±  5  meV 

Egx  -  hzzDBE 

DBE  localization  energy 

116meV  DBE  =  ZPL 

44  meV  DBE  =  1  —  P  replica 

Ed 

Donor  binding  energy 

25  ±  lOmeV 

Ea 

Acceptor  binding  energy 

175  ±25  meV 

255  ±  25  meV 

291  ±25  meV 

Epi 

MC  phonon — SBE 

52.5  meV/423  cm-1 

Ep2 

MC  phonon — DBE 

71.6meV/578  cm-1 
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FIG.  6.  (a)  Temperature  dependence  of  the  positions  of  the  FE-1P  [hz^FE-ip 
(T),  filled  triangles],  SBE-1P  [hi/BE_1P(T),  filled  squares],  and  DBE  (filled 
circles).  The  energy  separation  between  the  FE  and  SBE  lines  is  shown  as 
the  open  circles  and  is  shifted  up  by  3.27  eV  to  appear  on  the  same  scale. 
Subtracting  V^kT  from  this  (Eq.  (2))  gives  the  SBE  localization  energy  ABE 
(open  squares)  of  6.3  meV.  (b)  The  FE-1P  exciton  gap,  hi/pB-iFCT)  —  1/2kT,  is 
plotted  as  the  filled  circles,  and  is  fitted  with  both  Varshni  (dotted  line)  and 
Passler  (solid  line)  expressions.  This  fit  is  also  shown  shifted  up  by 
52.5  meV  (MC  phonon  energy)  to  represent  the  exciton  gap  Egx(T)  (dash-dot 
line)  and  by  an  additional  45  meV  (FE  binding  energy)  to  represent  the  indi¬ 
rect  band  gap  (dashed  line  at  the  top). 


open  circles.  These  are  shifted  up  by  3.27  eV  in  the  figure,  so 
as  to  appear  on  the  same  scale.  The  data  show  a  separation 
between  the  two  peaks  that  increases  slowly  with  tempera¬ 
ture,  as  expected  from  the  '/^kT  term  in  Eq.  (2).  ABE  is  deter¬ 
mined  by  subtracting  '/okT  from  this  data,  which  is  shown  as 
the  open  squares  and  results  in  a  temperature  independent 
BE  localization  energy  ABE  =  6.3  ±  1  meV. 

Following  the  FE  PL  up  to  room  temperature  provided 
data  that  could  be  fitted  over  almost  the  entire  temperature 
range  (25-294  K).  The  phonon-assisted  FE  gap,  fu/FE_1P(T) 
—  '/ikT,  is  plotted  as  the  filled  circles  in  Fig.  6(b)  and  was  fit¬ 
ted  with  the  ad  hoc  Varshni  equation,  s  Eee_  ip(T)  =  Efe-ip 
(0)  -  aT2/(T  +  b)  as  the  dotted  line  through  the  FE  data 
points,  with  EFE_1P(0)  =  3.372 eV,  a=-4x  10-4  eV/K  and 
b  =  — 640  K.  IBA  exhibits  negative  Varshni  coefficients, 
similar  to  some  other  wide  bandgap  semiconductors,  such  as 
6H-SiC  and  diamond.39  While  the  empirical  Varshni  equa¬ 
tion  has  been  the  expression  most  commonly  used  to  fit  the 
temperature  dependent  band  gap,  more  accurate  models  have 
been  proposed  that  have  the  advantage  of  depending  on  real 


physical  parameters.  A  model  developed  by  Passler40  has 
recently  been  employed  successfully  by  Grivickas  et  al 41  to 
fit  the  temperature  dependent  absorption  edge  of  4H-SiC. 
While  the  fonnal  Passler  expression  is  somewhat  compli¬ 
cated,  in  the  low-dispersion  limit  a  very  good  approxima¬ 
tion42  takes  the  form  E(T)  =  E(O)-(a0p/2){  [1  +  (2T/0p)p]1/p 
—  1 ) .  Here,  a  is  the  slope  of  E  vs.  T  at  the  high  temperature 
limit,  @p  is  the  average  phonon  temperature,  and  p  is  related 
to  the  degree  of  phonon  dispersion  A  by  p  =  (1  +  1/A“) /2.  The 
fit  of  the  FE  data  in  Fig.  6(b)  using  this  model,  which  was 
essentially  identical  to  that  using  the  full  expression,  is  given 
as  the  solid  line  through  the  FE  data  points,  with  Efe_1p  (0) 
=  3.372  eV,  oc  =  2x  10“3  eV/K,  0p  =  9OOK,  and  A  =  0.40 
(or  p  =  2.7).  At  low  temperature  and  room  temperature,  the 
fit  to  the  data  gave  Efe_1p  (0)  =  3.372  eV  and  Efe_1p  (294  K) 
=  3.275  eV. 

Except  for  Eg(0),  the  Passler  parameters  for  the  indirect 
band  gap  also  apply  to  those  determined  for  the  FE-1P  data. 
By  shifting  this  fit  of  the  FE  data  up  by  one  MC  phonon 
energy  (52.5  meV)  (see,  e.g.,  Eq.  (1)),  the  temperature  de¬ 
pendence  of  the  exciton  band  gap  Egx(T)  is  obtained  and  is 
given  by  the  dotted-dashed  line  in  the  figure.  Adding  to  this, 
the  FE  binding  energy  of  45  ±  5  meV  gives  the  temperature 
dependence  of  the  indirect  band  gap,  Eg(T),  indicated  in 
Table  II  and  shown  as  the  dashed  curve  at  the  top  of  the  fig¬ 
ure.  From  this  fit,  the  indirect  energy  gap  of  IBA  at  T  =  0 
and  at  room  temperature  is  obtained: 

E„(0)  =  3.470±0.005eV, 

(6) 

Eg (294  K)  =  3.373±0.01  eV. 

Since  the  fitted  FE  binding  energy  involves  a  5  meV  uncer¬ 
tainty  and  contributes  additively  to  the  low  temperature  band 
gap,  this  uncertainty  also  applies  to  Eg(0).  However,  this  cor¬ 
responds  to  a  rigid  energy  shift  of  Eg(T)  and  does  not  signifi¬ 
cantly  affect  the  temperature  dependence.  Obtaining  the 
room  temperature  band  gap  required  a  lit  of  the  data  to  the 
Passler  expression  and  a  total  uncertainty  in  Eg(294K)  of 
<10meV  is  estimated  at  room  temperature. 

The  «3.37  ±  0.01  eV  room-temperature  band  gap  appears 
in  between  the  two  previous  reports  of  3.20 eV  and  3.47  eV. 
The  fonner  was  determined  primarily  from  a  band  structure 
calculation,  as  the  ellipsometry  measurements  near  the  indi¬ 
rect  band  edge  were  dominated  by  impurity  absorption.21  This 
value  of  3.20 eV  for  the  indirect  band  gap  is  inconsistent  with 
the  present  PL  data,  however,  as  the  entire  294  K  near- 
bandedge  PL  spectrum  in  Fig.  4  lies  well  above  this  energy. 
The  100  me V  discrepancy  between  the  PL  and  the  optical 
absorption  results  is  rather  large.  To  investigate  this  further  the 
published  data  of  Slack  et  al24  have  been  digitized.  This  data 
were  originally  presented  in  the  fonn  of  absorption  coefficient 
vs.  photon  energy:  a  vs.  In/.  However,  indirect  absorption  is 
usually  plotted  as  a^2  vs.  \w  to  reflect  a  linear  dependence.  The 
original  data  of  Ref.  24  has  been  re-plotted  in  this  format  as 
the  heavy  solid  line  in  Fig.  7  and  the  positions  of  the  three  pro¬ 
posed  indirect  gap  energies  are  indicated.  The  dashed  line  is 
our  linear  extrapolation  of  the  absorption  data  that  leads  to  an 
intercept  with  the  energy  axis  at  3.47  eV,  similar  to  what  was 
obtained  by  Slack  et  al.  There  is,  however,  a  low  energy  tail, 
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hv  (eV) 

FIG.  7.  Digitized  optical  absorption  data  of  Ref.  24  (solid  line),  plotted  as 
a1'2  vs.  hi/.  The  dashed  line  represents  a  linear  extrapolation  with  a  threshold 
at  3.47 eV,  and  the  three  indirect  band  gaps  under  consideration  are  indi¬ 
cated  by  the  arrows.  An  extrapolation  of  the  tail  region  is  given  by  the  dotted 
line. 


which  might  be  interpreted  as  either  impurity  absorption,  indi¬ 
rect  absorption,  or  a  combination  of  the  two.  A  linear  extrapo¬ 
lation  of  the  tail  region  (dotted  line),  assuming  indirect 
absorption,  leads  to  an  intercept  at  about  3.33  eV,  well  below 
the  3.37  eV  band  gap  determined  here.  Consequently,  we 
believe  this  absorption  tail  to  be  primarily  due  to  impurity 
absorption,  similar  to  that  encountered  by  Bakalova  et  a/.,21 
which  masks  the  indirect  absorption  at  hv  >  3.37  eV.  The 
main  absorption  threshold  fitted  by  Slack  et  al.  occurs  at 
3.47  eV,  which  is  in  excellent  agreement  with  the  lowest  direct 
absorption  at  3.46 eV  reported  by  Bakalova  et  al21  While  this 
absorption  does  appear  to  exhibit  a  linear  tx^2  vs.  hv  depend¬ 
ence  above  about  3.6  eV,  the  available  evidence  suggests  that 
the  original  determination  of  the  indirect  absorption  edge  was 
probably  fitted  to  a  direct  transition.  In  their  original  paper, 
Slack  et  al.  stated  that  impurity  absorption  was  assumed  negli¬ 
gible.  However,  the  existence  of  a  tail  region  (Fig.  7)  should 
have  clearly  indicated  either  impurity  absorption  or  a  much 
lower  indirect  band  edge. 

Note  that  if  the  main  SBE  and  associated  FE  PL  bands 
were  interpreted  as  ZPLs  instead  of  lP-replicas,  so  that 
a  MC  phonon  energy  was  not  added  to  the  measured  FE 
line  positions  in  determining  Eg(T),  then  a  room  temperature 
band  gap  of  «3.32eV  would  have  been  obtained.  This  would 
be  in  excellent  agreement  with  the  extrapolated  tail  region  of 
the  Slack  et  al.  absorption  data,  and  the  absorption  tail  would 
consequently  be  assigned  as  the  indirect  absorption  edge. 
This  would  lead  to  a  «3.32eV  room  temperature  indirect 
band  gap  and  would  provide  a  consistent  explanation  of  both 
the  old  absorption  data  and  the  new  PL  data.  While  this 
would  make  for  a  very  satisfying  picture,  it  is  difficult  to 
ignore  the  correspondence  between  the  main  SBE  line  at 
3.366  eV  and  the  weak  3.419  eV  feature  in  Fig.  2(a),  which 
requires  that  the  dominant  BE  feature  be  assigned  as  a  1-P 
replica,  placing  the  room  temperature  band  edge  at  «3.37  eV. 

F.  Donor-acceptor  pair  spectra 

The  DAP  region,  ~3.0-3.3  eV,  is  dominated  by  three 
broad  peaks  centered  near  3.094,  3.13,  and  3.21  eV.  The  DAP 
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Energy  (eV) 

FIG.  8.  (a)  Dependence  of  the  DAP  spectrum  of  sample  A  on  excitation 
power.  The  peak  intensities  are  all  normalized  to  the  same  height,  (b)  The 
high  energy  edge  of  the  DAP  spectrum  shown  at  low  and  high  laser  power. 
The  sharp  and  reproducible  features  in  the  high  power  spectrum  reflect 
recombination  from  close  DAPs  at  discrete  donor-acceptor  separations. 


nature  of  the  emission  becomes  apparent  from  the  depend¬ 
ence  of  the  spectrum  on  excitation  power.  In  Fig.  8(a),  we 
compare  three  spectra  (with  all  peak  intensities  normalized  to 
the  same  amplitude)  obtained  using  relative  laser  powers  of 
1.0,  0.1,  and  0.01.  With  decreasing  power,  the  high  energy 
edge  of  the  two  higher  energy  lines  shifts  to  lower  energy — 
one  of  the  hallmarks  of  DAP  recombination.43  The  apparent 
peak  positions  are  affected  by  both  the  varying  high  energy 
edge  and  the  varying  sloped  backgrounds  that  are  due  to  the 
lower  energy  DAPs.  Consequently,  the  shifts  in  position  of 
some  of  the  peaks  do  not  follow  the  shift  in  the  high  energy 
edge  in  a  linear  fashion.  As  we  cannot  see  the  high  energy 
edge  for  the  lowest  3.094 eV  band,  we  can  only  assume  by  its 
proximity  to  the  other  lines  that  it  is  also  of  DAP  origin.  As 
the  three  lines  were  observed  to  exhibit  different  dependences 
upon  temperature  (see  Fig.  9),  they  do  not  appear  to  be  differ¬ 
ent  phonon  replicas  of  a  single  pair  band,  but  instead  repre¬ 
sent  independent  impurity  pairs.  As  there  has  not  been  a 
previous  report  of  a  donor  in  IBA,  it  is  most  likely  that  the 
three  lines  represent  three  different  acceptors.  Also,  since 
IBA  is  an  indirect  gap  material,  the  DAP  emission  could  also 
involve  a  momentum-conserving  phonon.  However,  no  asso¬ 
ciated  ZPL  was  identified. 
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FIG.  9.  (a)  Temperature  dependence  of  the  DAP  spectrum  between  10  K 
and  110K.  Peak  intensities  are  all  normalized  to  the  same  height,  corre¬ 
sponding  to  the  intensity  of  the  deepest  DAP  feature  at  3.094  eV,  except  for 
the  110K  data.  The  intensities  of  the  higher  energy  lines,  corresponding  to 
shallower  impurities,  are  seen  to  decrease  with  increasing  temperature  faster 
than  the  deeper  bands,  (b)  Temperature  dependence  of  the  integrated  inten¬ 
sity  of  the  3.21  eV  band.  The  solid  line  is  a  fit  of  the  data  to  Eq.  (8). 


The  high-energy  edge  of  the  DAP  region  of  sample  A  is 
shown  in  Fig.  8(b)  for  two  power  levels  differing  by  an  order 
of  magnitude.  A  series  of  sharp,  reproducible  emission  bands 
extending  from  the  broad  distant-pair  band  («3.21eV)  to 
higher  energy  are  apparent  at  high  excitation  and  are  consid¬ 
erably  reduced  at  lower  power.  These  represent  recombina¬ 
tion  between  close  donor-acceptor  pairs  at  discrete  DAP 
separations,43  similar  to  those  observed  in  BP,44  SiC,45  and 
GaP.46  The  observation  of  a  close-pair  spectrum  further  con¬ 
firms  the  DAP  origin  of  the  3. 0-3. 3  eV  region.  While  we  can¬ 
not  yet  identify  the  donor  and  acceptor  species  responsible 
for  the  observed  PL  emission,  the  existence  of  a  DAP  band 
provides  evidence  of  a  donor  in  IBA,  which  is  of  consider¬ 
able  interest,  as  n-type  material  has  not  yet  been  reported. 

In  /1-rhombohedral  boron,  Ni  as  well  as  several  other 
transition  metals  were  shown  to  be  effective  as  n-type  dop¬ 
ants47'48  that  incorporate  on  an  interstitial  site  of  the  boron 
icosahedron.48  Ni  also  acts  as  an  n-type  dopant  in  boron  car¬ 
bide49  and  in  boron-carbon  alloys.50  Given  the  high  concen¬ 
tration  of  Ni  in  these  IBA  samples  (Table  I)  and  the  fact  that 
the  structure  of  IBA  shares  the  boron  icosahedron  with  these 


compounds,  it  is  likely  that  the  donor  in  question  is  an  inter¬ 
stitial  Ni  impurity. 

The  possibility  of  a  Ni  donor  in  IBA  is  discussed  in 
more  detail  below.  However,  since  a  Ni-related  DBE  was  al¬ 
ready  suggested  in  an  earlier  section,  some  further  comment 
seems  appropriate.  For  both  the  shallow  donor  and  the  DBE, 
it  should  be  made  clear  that  no  evidence  is  presented  that 
definitely  identifies  Ni  with  either  of  these  two  centers.  There 
is  sufficient  circumstantial  evidence  presented,  however,  to 
justify  serious  consideration  of  the  role  of  Ni  in  both  of  these 
cases.  The  shallow  donor  and  the  DBE  would  correspond  to 
Ni  incorporation  at  different  sites,  which  would  certainly  not 
be  surprising,  given  the  large  Ni  concentration  in  these 
samples. 

Since  electrical  measurements  have  not  been  carried  out 
in  the  present  samples,  and  with  the  high  concentration  of  Ni 
impurities,  we  cannot  be  sure  if  the  samples  are  p-type,  as 
they  usually  are,  or  n-type.  We  also  have  no  direct  informa¬ 
tion  about  the  compensation  level  in  these  samples.  How¬ 
ever,  the  observation  of  relatively  sharp  exciton  structure  at 
low  temperature  is  not  consistent  with  a  high  concentration 
of  free  carriers.  This  suggests  that  the  material  is  highly 
compensated  and/or  that  only  a  small  fraction  of  the  Ni 
impurities  incorporate  as  interstitial  donors  in  the  boron  ico- 
sahedra.  Also,  since  the  PL  measurements  were  restricted  to 
the  high  quality  region  of  each  sample,  the  Ni  concentrations 
in  these  regions  may  be  much  lower  than  what  was  indicated 
in  Table  I.  Evidence  of  significant  compensation  is  observed, 
however,  in  the  SBE  lineshape  in  Fig.  2(a),  which  exhibits  a 
clear  asymmetric  tail  at  lower  energies.  This  BE  lineshape 
has  been  associated  both  theoretically51  and  experimen¬ 
tally52  with  the  Stark  effect  due  to  ionized  impurities.  The 
lineshape  in  Fig.  2(a)  could  be  fitted  to  the  model  in  Ref.  51, 
which  provided  a  rough  estimate  of  the  ionized  impurity  con¬ 
centration,  w2x  1018cm" 3.  However,  in  order  to  confirm 
the  role  of  Ni  in  the  transport  properties  of  IBA,  a  combina¬ 
tion  of  electrical  and  optical  measurements  on  samples  with 
carefully  controlled  Ni  concentrations  will  be  required. 

To  obtain  an  estimate  of  the  donor  and  acceptor  binding 
energies,  the  DAP  spectrum  was  studied  as  a  function  of 
temperature.  In  general,  the  measured  photon  energy  h//|)AP 
emitted  by  an  indirect  band  gap  material  through  a  DAP 
transition  between  neutral  donors  and  acceptors  separated  by 
a  distance  r  is  given  by 

ht'DAp(r)  =  Eg-(ED  +  Ea)  +  e2/ sr  -  EP,  (7) 

where  Eg,  ED,  and  EA  are  the  band  gap  and  donor  and 
acceptor  binding  energies,  respectively,  e2/sr  is  the  Coulomb 
interaction  between  the  two  ionized  impurities  of  the  final 
state,  with  e  the  electronic  charge  and  e  the  static  dielectric 
constant.  EP  is  the  energy  of  the  MC  phonon.  Fig.  9(a)  shows 
DAP  spectra  taken  between  10  K  and  110K,  with  the  inten¬ 
sities  (monotonically  decreasing  with  temperature)  all  nor¬ 
malized  to  the  same  peak  height.  It  is  apparent,  as  noted 
earlier,  that  the  3.21  eV,  3.13  eV,  and  3.094  eV  peaks  exhibit 
different  dependences  on  temperature,  with  the  highest 
energy  emissions  decreasing  more  rapidly  with  temperature 
than  the  deeper  bands.  This  is  exactly  what  is  expected  from 
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the  temperature  dependence  of  DAP  PL  involving  more  than 
one  acceptor  (or  donor):  The  highest  energy  emission  corre¬ 
sponds  to  the  shallowest  acceptor  (donor),  which  thermalizes 
first,  before  the  deeper  impurities. 

At  higher  temperatures  (>90K),  a  broad  band  centered 
near  3.0  eV  and  with  intensity  decreasing  much  more  slowly 
with  temperature  than  the  DAP  band  emerged  from  the  low- 
energy  tail  of  the  DAP  and  had  a  significant  effect  on  the 
low  energy  portion  of  the  DAP  spectrum  at  these  tempera¬ 
tures.  This  band  totally  dominates  the  1 10  K  spectrum  in  Fig. 
9(a).  While  its  origin  is  unidentified  at  present,  the  large  line 
width  suggests  recombination  at  a  defect  center. 

To  extract  the  temperature  dependence  of  the  DAP 
intensities,  it  is  necessary  to  separate  the  contributions  from 
the  three  individual  peaks  (as  well  as  the  3.0  eV  band). 
Because  the  two  lowest  energy  peaks  in  the  DAP  spectra  are 
relatively  close  together,  it  was  not  possible  to  reliably 
deconvolute  these  bands.  However,  the  3.21  eV  peak  is  well 
separated  from  the  others  and  could  be  successfully  sepa¬ 
rated  from  the  background  of  the  other  lines,  which  allowed 
its  integrated  intensity  to  be  determined  as  a  function  of  tem¬ 
perature.  This  dependence  is  shown  as  the  filled  circles  in 
Fig.  9(b).  The  data  are  fitted  with  two  ionization  paths,  corre¬ 
sponding  to  the  sequential  thermalization  of  the  donor  and 
the  acceptor: 

Idap(T)  oc  t  +  Ae-ED/kT  +  Be-EA/kT  ■  (8) 

The  best  fit  to  the  data  in  Fig.  9(b),  assuming  that  Ea  >  ED, 
was  obtained  with  EA=  175  ±  25meV,  ED  =  25  ±  lOmeV. 
Using  the  measured  separations  between  the  3.21  eV  band 
and  the  other  two  deeper  DAP  lines,  we  estimate  acceptor 
binding  energies  of  255  meV  and  291  me V  for  the  two 
deeper  acceptors. 

Recent  Hall  effect  measurements15  on  heteroepitaxial 
layers  of  IBA  on  semi-insulating  6H-SiC  substrates  found 
the  as-grown  layers  to  be  p-type.  Three  shallow  acceptor 
activation  energies  were  reported  at  150,  260,  and  330  meV. 
This  agrees  well  with  the  acceptor  binding  energies  of  175, 
255,  and  291  me V  obtained  from  the  DAP  PL  spectra  and 
further  supports  a  DAP  band  derived  from  three  distinct 
acceptors  and  a  single  donor.  That  the  donor  is  the  shallower 
of  the  two  impurities  is  also  consistent  with  the  electron 
effective  mass  calculated  by  Bakalova  el  a/.21  (me*  =  0.29), 
which  was  considerably  less  than  the  corresponding  hole 
masses.  It  is  also  noteworthy  that  the  fitted  donor  binding 
energy  of  25meV  is  substantially  lower  than  that  expected 
from  an  effective  mass  donor:  ED  =  me*RH/e2  =  64  meV, 
where  me*  and  RH  represent  the  electron  effective  mass21 
and  the  hydrogen  Rydberg,  13.6  eV,  respectively,  and 
e  =  7.84  (Ref.  21).  Such  a  small  binding  energy  is  consistent 
with  an  interstitial  donor  impurity,  such  as  the  Ni  impurity 
suggested  above.  Interstitial  shallow  impurities  with  binding 
energies  considerably  smaller  than  that  predicted  by  effec¬ 
tive  mass  theory  have  been  reported  in  the  literature.  Some 
examples  are  Li  and  Na  interstitials  in  ZnSe,5’  Cl  interstitials 
in  CuGaS2,54  and  Zn  interstitials  in  ZnQ.  In  the  case  of  the 
latter,  the  deviation  from  effective  mass  theory  was  associ¬ 


ated  with  the  non-spherical  charge  distribution  around  the 
impurity  site. 

Thus,  the  temperature  dependence  of  the  intensity  of  the 
shallowest  DAP  band  at  3.21  eV  corresponds  to  shallow 
impurities  with  ED  +  EA  «  200  ±  30  meV.  We  can  make  an 
independent  estimate  of  ED  +  EA  by  employing  Eq.  (7), 
using  Eg  =  3.47eV,  h^DAP=  3.21  eV,  and  EP=52.5meV. 
Dean4’  has  given  an  estimate  of  the  DAP  separation  corre¬ 
sponding  to  the  peak  of  the  distant  DAP  band: 
rpk  =  ( '/27rN ) _  1  ''  L  where  N  is  the  concentration  of  the  domi¬ 
nant  impurity.  For  N  in  the  range  1 0 1 6—  1 0 1 8  cm-3,  the  corre¬ 
sponding  Coulomb  energy  is  in  the  range  (5-21)  meV. 
Equation  (7)  then  gives:  (ED  +  EA)  =  (213-229)  meV  for 
Ep  =  52.5meV  and  (ED  +  EA)  =  (265-281)  meV  for  EP  =  0 
(i.e.,  if  no  MC  phonon  is  involved).  The  experimental  result 
of  (Ed  +  Ea)  =  200  ±  30  meV  =  (170-230)  meV  is  consist¬ 
ent  with  a  DAP  emission  that  is  phonon-assisted ,  employing 
a  52.5  meV  phonon.  While  no  ZPLs  were  identified  for  the 
three  DAP  bands,  the  ZPL  positions  would  be  projected  to 
be  at  3.146,  3.182,  and  3.262eV. 

Within  the  effective  mass  approximation,  the  acceptor 
binding  energy  is  related  to  the  hole  effective  mass  through 
EA  =  mh*RH/e2.  With  the  measured  acceptor  binding  energy 
for  the  3.21  eV  band  in  the  range  EA=  (150-200)  meV,  we 
find  the  corresponding  mh*  in  the  range  (0.68-0.90).  This  is 
in  stark  contrast  to  the  huge  heavy-hole  effective  mass  of 
mhh*  =  359  calculated  for  the  I  -point  in  Ref.  21.  Combining 
mh*  =  0.68-0.90  with  me*  =  0.29,  we  can  obtain  an  estimate 
of  the  reduced  mass  mr*  necessary  to  calculate  the  FE  bind¬ 
ing  energy,  l/mr*  =  l/me*  +  l/mh*,  with  the  result  that 
mr*  =  (0.203-0.219).  The  corresponding  FE  binding  energy, 
AFE  =  mr*  RH/e2,  is  then  in  the  range  (45-48.5)  meV,  which 
is  in  much  better  agreement  with  the  measured  FE  binding 
energy  than  the  66meV  calculated  in  Ref.  36  assuming  an 
infinite  hole  mass. 

The  experimental  data  indicate  that  the  DAP  spectrum 
involves  relatively  shallow  acceptors  and  that  the  FE  binding 
energy  estimated  above  from  the  shallow  acceptor  binding 
energy  is  in  good  agreement  with  experiment.  On  the  other 
hand,  the  very  large  heavy-hole  effective  mass  determined 
from  the  band  structure  calculations  simply  reflects  a  very 
flat  valence  band.  This  was  also  a  feature  of  several  earlier 
band  structure  calculations  for  IBA  (Refs.  20-22)  and  is  pre¬ 
sumably  accurate.  However,  such  a  large  hole  mass  cannot 
support  a  shallow  acceptor.  This  suggests  that  the  observed 
shallow  acceptors  might  be  associated  with  the  light-hole 
band.  Bakalova  et  a/.21  calculated  a  light-hole  mass  of 
m]h*  =  0.46.  This  corresponds  to  an  effective  mass  acceptor 
binding  energy  of  102  me V.  Given  that  mlh*  has  not  been 
measured  directly,  and  that  the  fitted  binding  energy  of  the 
shallowest  acceptor,  150-200  meV,  undoubtedly  contains 
some  central  cell  contribution,  this  could  be  considered  as 
reasonable  agreement.  Using  the  calculated  light-hole  mass 
gives  a  reduced  mass  for  the  FE  of  mr*  =  0.178  and  a  FE 
binding  energy  of  AFE  =  39.4  me  V,  which  is  also  in  reasona¬ 
ble  agreement  with  experiment.  Consequently,  we  conclude 
that,  due  to  the  very  large  heavy-hole  mass,  the  observed 
shallow  acceptors  and  the  FE  are  probably  associated  with 
the  light-hole  band. 
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IV.  SUMMARY 

PL  spectroscopy  was  carried  out  on  solution-grown  IBA 
crystals  in  the  range  1.75-3.5  eV  and  over  the  temperature 
range  8K-294K  in  order  to  measure  the  indirect  band  gap 
and  to  investigate  the  shallow  impurities  in  the  material.  The 
spectra  were  dominated  by  recombination  due  to  deep 
defects,  shallow  donor-acceptor  pairs,  and  excitons  bound  to 
both  deep  and  shallow  centers.  Analysis  of  the  temperature 
dependence  of  the  PL  intensities  and  line  positions  provided 
several  material  parameters,  including  the  temperature  de¬ 
pendence  of  the  energy  gap  and  the  binding  energies  of  the 
excitons  and  shallow  donors  and  acceptors.  These  are  sum¬ 
marized  in  Table  II. 

The  dominant  SBE  emission  was  shown  to  involve  a 
phonon-assisted  transition,  while  the  strongest  DBE  line 
appeared  to  be  a  ZPL,  as  no  ZPL  at  higher  energy  could  be 
observed.  FE  recombination  was  observed  for  temperatures 
>25  K  and  became  the  dominant  feature  in  the  spectrum  by 
120  K.  The  temperature  dependence  of  the  PL  intensity  of 
the  FE  provided  the  FE  binding  energy  (45  meV),  which  was 
a  necessary  parameter  to  determine  the  indirect  band  gap. 

The  shift  of  the  energy  position  of  the  phonon-assisted 
FE  line  with  temperature  reflects  the  temperature  depend¬ 
ence  of  the  indirect  band  gap  and  was  fitted  with  both  Var- 
shni  and  Passler  expressions.  The  resulting  indirect  band  gap 
at  low-  and  room-temperatures  was  determined  as  3.470  eV 
and  3.373  eV,  respectively.  This  room  temperature  band  gap 
is  in  contrast  to  two  previously  published  values  for  the  indi¬ 
rect  gap:  3.20  eV  (band  structure  calculation)  and  3.47  eV 
(optical  absorption).  The  former  was  found  inconsistent  with 
the  observed  room  temperature  PL  spectra.  A  re-examination 
of  the  original  optical  absorption  spectra,  coupled  with 
recent  ellipsometry  measurements,  suggests  that  the  optical 
absorption  data  appear  to  have  been  fitted  to  the  lowest 
direct  transition  rather  than  the  indirect  edge,  which  may 
have  been  masked  by  impurity  absorption. 

A  strong  deep  bound  exciton  emission  was  observed  at 
3.309  eV.  The  temperature  dependence  of  the  DBE  spectrum 
reflected  a  peak  position  that  was  independent  of  temperature 
up  to  >100K.  This  behavior  is  reminiscent  of  earlier  reports 
of  DBEs  bound  to  transition  metals  and  suggests  the  possible 
involvement  of  a  Ni-related  site  in  the  origin  of  the  DBE. 

DAP  spectra  were  identified  through  their  dependence 
on  excitation  power,  temperature  and  by  the  appearance  of 
discrete  DAP  emission  from  close  pairs.  Three  distant  DAP 
bands  were  observed  and  were  associated  with  three  differ¬ 
ent  acceptors.  The  observation  of  DAP  emission  in  IBA  indi¬ 
cates  the  presence  of  a  donor.  This  is  a  significant  finding,  as 
n-type  material  has  not  yet  been  reported.  While  the  identity 
of  the  donor  impurity  remains  unknown,  it  is  suggested  that 
the  donor  may  result  from  a  fraction  of  the  large  Ni  concen¬ 
tration,  incorporated  on  interstitial  sites  of  the  boron  icosahe¬ 
dron.  Ni  doping  has  produced  n-type  material  (via  a  Ni 
interstitial)  in  other  compounds  incorporating  boron  icosahe- 
dra,  and  the  small  observed  donor  binding  energy  observed 
in  IBA  is  consistent  with  an  interstitial  donor.  The  tempera¬ 
ture  dependence  of  the  intensity  of  the  shallowest  DAP  band 
provided  estimates  of  the  donor  and  acceptor  binding  ener¬ 


gies  (see  Table  II).  However,  the  very  large  heavy-hole 
effective  mass  reported  from  recent  band  structure  calcula¬ 
tions  cannot  support  such  shallow  acceptors.  The  results  are, 
therefore,  interpreted  as  due  to  acceptors  associated  with  the 
light-hole  band. 
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